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We analyze two-photon double ionization of helium in both the nonsequential (hw < I2 ~ 54.4 eV) 
and sequential {fwj > I2) regime. We show that the energy spacing A_B = E\ — E2 between the 
two emitted electrons provides the key parameter that controls both the energy and the angular 
distribution and reveals the universal features present in both the nonsequential and sequential 
regime. This universality, i.e., independence of hjj, is a manifestation of the continuity across the 
threshold for sequential double ionization. For all photon energies, the energy distribution can be 
described by a universal shape function that contains only the spectral and temporal information 
entering second-order time-dependent perturbation theory. Angular correlations and distributions 
are found to be more sensitive to the value of hio. In particular, shake-up interferences have a 
large effect on the angular distribution. Energy spectra, angular distributions parameterized by 
the anisotropy parameters , and total cross sections presented in this paper are obtained by fully 
correlated time-dependent ab initio calculations. 

PACS numbers: 32.80.Rm, 32.80.Fb, 42.50.Hz 



I. INTRODUCTION 

Two-photon double-ionization (TPDI) of helium is a 
prototype process for the study of electron correlation. 
The advent of ultrashort and intense light sources with 
sufficient photon flux [HI] offers the opportunity to in- 
vestigate TPDI over a wide range of energies, from direct 
(nonsequential) to sequential ionization. For photon en- 
ergies below I2 ~ 54.4 eV, the second ionization potential 
of helium, one photon does not carry sufficient energy to 
ionize the residual He"*" ion in its ground state. After 
the first ionization event below this threshold, the sec- 
ond photoabsorption event has to occur within a short 
time interval such that the outgoing electrons can ex- 
change energy for the double ionization process to take 
place. 

However, above this threshold, sequential double ion- 
ization (SDI) becomes possible: the first photon singly 
ionizes the helium atom and the second photon ionizes 
the remaining He"'' ion, each of which constitutes a sepa- 
rate on-shell process. The time interval elapsed between 
the two photoabsorption events and, thus, between the 
electron emission events, can be, in principle, arbitrarily 
long. This sequential process can be qualitatively well 
described within an independent particle picture while 
quantitative details are infiuenced by electron-electron 
interactions [12-23]. 

By contrast, electron correlation is a conditio sine 
qua non for the direct nonsequential double ionization 
(NSDI) process to occur. Therefore, much effort has 



* renate.pazourek@tuwien.ac.at 



been spent on investigations of electron correlation in 
the nonsequential regime. So far, partially integrated 
quantities and total cross sections have been measured 
[5, 24-27]. Several theoretical studies of fully differential 
cross sections have been presented [12-14, 28-30]. 
However, even on the level of total cross sections the 
different theoretical approaches lead to differences of 
more than one order of magnitude ]12, 14-17, 30-32]. 
The reasons for these discrepancies are still under 
debate, and the experiments could, up to now, not reach 
the needed accuracy to resolve them. The extraction 
of cross sections close to the threshold for sequential 
ionization has remained a challenging problem. 

The focus of the present communication is on the close 
similarity and underlying common features of the NSDI 
below the threshold I2 and the SDI above the threshold 
l2- While the continuity across thresholds is appreciated 
as a consequence of the analyticity of the S-matrix and 
is frequently involved in the determination of threshold 
exponents (Wigner and Wannier exponents [33, 34]), its 
implication for energy- and angular distributions in the 
TPDI process has, so far, not been explored. 

We show that energy distributions and, to a consid- 
erable extent, also angular distributions of TPDI dis- 
play universal features present both above and below the 
threshold for SDI. These universal features become obvi- 
ous when observables are analyzed in terms of the energy 
spacing (or asymmetry of energy sharing) /S.E — Ei — E2 
of the two outgoing electrons. The significance of this pa- 
rameter controlling the double ionization irrespective of 
the value of huj can be understood within the framework 
of second-order time-dependent perturbation theory. We 
will discuss the scaling of the resulting shape function 
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with pulse duration T. 

Application to angular distributions characterized by 
anisotropy parameters fij shows that /S.E still plays a 
significant role even though the j3j refer to reduced one- 
electron variables. The latter are of experimental rel- 
evance, as a coincidence measurement is not required, 
at least for electron energies where the spectrum is not 
dominated by single ionization. We compare the single- 
electron angular distribution of TPDI with some results 
from earlier calculations [18, 19, 35], and find significant 
differences. 

Unless otherwise stated, atomic units will be used 
throughout the text. 



II. METHOD 

In our computational approach (see [30, 36, 37] for a 
more detailed description) we solve the time-dependent 
Schrodinger equation in its full dimensionality, includ- 
ing all interparticle interactions. We employ a time- 
dependent close-coupling scheme ]20, 28, 29, 38] where 
the angular part of the wave function is expanded in cou- 
pled spherical harmonics. In order to reach convergence 
in the angular coordinates, we use single-electron angular 
momenta up to values of /i^max = ^2. max = 10. The high- 
est total angular momentum included in the time propa- 
gation is typically imax — 2, which is sufficient since only 
the two-photon channels L ~ Q and L = 2 play a signifi- 
cant role. We checked explicitly that using imax — 3 does 
not change the results presented providing clear evidence 
that lowest-order perturbation theory in the photon field 
provides the dominant contribution. For the discretiza- 
tion of the two radial variables, we employ a finite ele- 
ment discrete variable representation (FEDVR) ]39-42]. 
The sparse structure of the resulting matrices allows for 
efficient parallelization, which is crucial to obtain results 
in the long-pulse limit (up to more than 20 femtoseconds 
propagation time). Radial boxes with an extension of 
up to Tinax = 2000 a.u. Containing FEDVR elements with 
lengths of 4 — 4.4 a.u. and of order 11 are used for the pre- 
sented results. For the temporal propagation of the wave 
function, we employ the short iterative Lanczos method 
[43-45] with adaptive time-step control. 

The laser field is assumed to be linearly polarized and 
treated in dipole approximation. The interaction oper- 
ator is implemented in both length and velocity gauge, 
such that gauge independence can be explicitly verified. 
We choose the temporal shape of the vector potential to 
be of the form 

A{t)^zAosm^ (^^^siniujt) (1) 

for < t < 2T. The FWHM of the sin^ envelope function 
has the duration T. A peak intensity of Iq — 10^^ W/cm^ 
ensures that ground state depletion and three- or higher- 
order photon effects are negligible. 



The asymptotic momentum distribution is obtained by 
projecting the wave packet onto products of Coulomb 
continuum states. These independent-particle Coulomb 
wave functions are not solutions of the full Hamiltonian. 
However, as we have previously demonstrated, projection 
errors can be reduced and controlled to the one-percent 
level by delaying the time of projection until the two elec- 
trons are sufficiently far apart from each other ]30]. All 
results were tested for numerical convergence and gauge 
independence. 



III. SHAPE FUNCTION FOR TWO-PHOTON 
DOUBLE IONIZATION 

The point of departure of our analysis of common fea- 
tures of TPDI above and below the sequential thresh- 
old observed in the numerical calculations is the spec- 
tral shape function within second-order time dependent 
perturbation theory. The applicability of the latter is 
indicated by the negligibly small contributions of three- 
photon (or higher-order) processes. Similar approaches 
have been employed previously, see ]12, 16] and references 
therein. To second order the amplitude of the transition 
driven by the xuv pulse is 

*/ tl 

= - ^ J<itil dt2e^^^"*^e'^-*^ X 

n to to 

X (/|V(ti)|n)(n|V(i2)N) (2) 

with Efn ~ Ef — En and Eni = En — The transition 
probability for TPDI is P^^ {Ei,ni, £2,^2) = 1*^/^, 
where \i) is the ground state and |/) = \EiQi, E2i^2)- 
Insertion of the interaction operator in velocity gauge, 
V(i) = (pz,i +i>z,2)A{t) = fiA{t), leads to a factorization 
of each term in the sum over intermediate states into a 
spectral function Q that just depends on the energies of 
the involved states and the temporal shape of the inter- 
action potential, and a time independent matrix element 
depending on two dipole operators, 

^-^{f\f^\n){nm)giEfn,E„,,[Am (3a) 

n 

with 

t/ fi 

g{Efn,Er,.„[Am= Jdti Jdt2e^''f-'^e'^-^'^A{t,)A{t2). 

to to 

(3b) 

Employing the rotating wave approximation and using 
the temporal shape of the pulse (Eq. 1), Eq. 3b becomes 
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Al r f 4a 
I dh I dt2FiAEf„,h,T)F{AEm,t2,T) 



with 

F{6,t,T) ^e^'' sin' (^^y (4b) 

where AEfn = Efn — w and Ai?„j = — uj. The latter 
variables denote the energy defects relative to the energy 
conserving on-shell transition in each of the two steps. 
The integral can be solved analytically, but is not shown 
here for brevity. 

Above the sequential threshold {hw > I2) and for 
infinitely long pulses T — ^ 00 the angle-integrated 
joint transition probability P^^ {Ei, E2) is governed by 
G{AEf„ = 0,AEni = 0,Ao,T), i.e., the on-shell part 
of the transition amplitude. For short pulses and/or 
Hll! < I2, the behavior of P^^ {Ei, E2) is controlled by 
the deviations AEfn 7^ and AEni 7^ 0. Eq. 3a can 
be further simplified by noting that among all possible 
intermediate states, the continuum states \Ep, Is) associ- 
ated with the ionic ground state and an outgoing p-wave 
strongly dominate. In other words, shake-up or shake- off 
processes only represent a few-percent correction in the 
SDI regime. 

We further make the approximation that the transi- 
tion amplitude from the intermediate to the final state is 
diagonal in the energy Ei of the free electron, i.e., that 
{Eini,E2n2\fi\E[p,ls) cx S{E[ - El). Consequently, a 
single term |no) = \Eip, Is) in the sum over intermediate 
states provides the leading contribution. 

Instead of using an independent-particle approxima- 
tion for the transition matrix elements (as in [12, 16]), 
we approximate them by a constant. This can be done if 
one is only interested in the shape of the final electron en- 
ergy distribution. The transition probability to the final 
energies {Ei,E2) in this approximation is given by 

P^'{Ei,E2) = C7|g3i„^(Ai?^„a), Ai?„(i),,T)+ 

^,i„2(Aii;^„(.,,A£;„(.,^,r)p, (5) 

with AE^ (1) = Ef - (El + + u). AE^ (2) = 

Ef - {E2 + E'-^"'^ + uj) and corresponding expressions for 
AE Ik) .. 

Accordingly, Eq. 5 can be rewritten as 

P^\Ei,E2) - Pi'{AE, Etot.i) , (6) 

with AE — El — E2- In the limit T — > 00, the total en- 
ergy of the final state -Etotai — E1+E2 = 2hu! -\- Eg is well 
determined by the photon energy and the ground state 
energy of helium, Eq sa — 79eV. For short pulses with 



finite Fourier width, the distribution of i?totai is broad- 
ened accordingly. The reduced probability density as a 
function of AE follows after integration over the Fourier 
width as 

P^'iAE) = P<f^(i?i,i?2)di?total (7) 

(with a factor 1/2 from the Jacobi determinant of the 
coordinate transformation). For long pulses, this projec- 
tion through the (i?i, E2) plane closely approximates the 
conventional one-electron energy distribution. Eq. 7 has 
a distinct advantage when comparing pulses with differ- 
ent photon energies: the sequential peaks at Huj — Ii and 
huj — 12 always show up at the same positions ±(/i — 12), 
irrespective of the photon energy. 

We refer to Eq. 5 as the shape function for double ion- 
ization. It contains both the spectral information of the 
time-independent fully correlated two-electron Hamilto- 
nian and the temporal structure of the pulse. It thus 
depends on the energy exchange between the two elec- 
trons included via ionization potentials that depend on 
the presence of the second electron, i.e., Ii 7^/2- How- 
ever, no additional information on electron correlations, 
in particular angular correlations, is included. 

In Fig. 1 we compare Pg-'{AE) (Eq. 7) with the nu- 
merically exact solution of the time-dependent Schro- 
dinger equation for a wide range of photon energies 
A2e\ < hio < 80 eV covering both the NSDI and SDI 
regimes with a pulse duration of T = 4.5 fs for all pulses. 
For better comparison of the results at different photon 
energies, we have normalized each curve to the value at 
AE = 0. 

The excellent agreement with Pg^{AE) for energy dif- 
ferences smaller than \AE\ Ri 45 eV irrespective of Itoj ex- 
hibits the universal features present both in the SDI and 
NSDI process (note that the deviations visible for larger 
\AE\ result from the neglect of intermediate shake-up 
states in the sum in Eq. 3a). Non- vanishing emission 
probabilities far away from the "sequential" peaks with 
AEfn — Ai?„i — signify energy sharing via electron- 
electron interaction irrespective of whether the sequential 
peaks are energetically accessible or not. Thus, also for 
pulses in the nominally sequential regime (huj > I2), a 
large part of the final-state space can only be accessed 
through energy sharing that is governed by electron- 
electron interaction and thus closely resembles the NSDI 
process. Only in the immediate vicinity of the on-shell 
peaks AEfn = ^Eni = 0, the width of which decreases 
as T~^, SDI and NSDI processes become distinct. 

IV. ANGULAR DISTRIBUTIONS 

We extend now the analysis of the DI probability as 
a function of AE to a second degree of freedom, the 
emission angle of one of the emitted electrons with re- 
spect to the laser polarization axis, 9. Conventionally the 
one-electron TPDI distribution P^^{E,9) is considered. 
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FIG. 1. (a) Singly differential energy distribution P^^ (AE) as a function of the energy difference AE — Ei ~ E2, normalized 
to the yield at AE = for photon energies between 42 eV and 80 eV. The pulses had a duration T = 4.5 fs and an intensity 
Jo — 10^^ W/cm^. The gray lines show the expected positions of the peaks for the sequential process for different intermediate 
states (i.e., with and without shake-up). Here, £„ is the excitation energy from the ionic n — 1 state to an excited state n 
{£2 ~ 40.8 eV). The black line shows the energy distribution based on second-order perturbation theory Pg^{AE). (b) and 
(c) show the two-electron energy distribution P^^ [Ei, E2) on the log scale for (b) hw — 48 eV and (c) hu — 80 eV. 




FIG. 2. (a) Anisotropy parameters j3{u},AE) and ^{lj,AE) (see Eq. 9) for photon energies between 42 eV and 80 eV as a 
function of the energy difference AE = Ei — E2 in the interval [—Etot, Etot], with Etot = 2ftcj + Eq, Eq ~ — 79eV. The pulse 
parameters are the same as in Fig. 1. The upper group of lines contains the /3 parameters, while the lower group contains the 
values for 7. Note that the angular distribution of the electron with Ei is shown, thus A_E > characterizes the faster electron 
and AE < the slower one. (b) Angular distribution at equal energy sharing P^^ {AE = 0, 9) for huj = 42 eV, 48 eV and 70 eV 
(from inside to outside, same color code as in (a)), normalized to a value of 1 for 6 = 0°. The laser polarization axis is indicated 
by the arrow. The thin black line shows a cos^(6l) distribution. 



which is obtained by integrating the full two-electron dis- 
tribution P'^'{Ei,E2,fli,n2) over the energy and angle 
of one electron. This corresponds to measuring electrons 
without coincidence requirements. Following our findings 
that the TPDI distribution depends primarily on AE, we 
therefore investigate in the following P^^ {AE,9i). Two 



features are worth noting: since AE' ss 2Ei — {Eci + 2hu>), 
the switch from Ei to AE corresponds to a coordinate 
shift in the limiting case of long pulses, i.e., narrow 
Fourier width of hio. Furthermore, since only the angle 
of one electron {9i) is observed, P^^ {AE,6i) is, unlike 
P^^ (AE), not symmetric under inversion AE — )■ —AE. 
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The one-electron distribution has cylindrical symmetry 
(i.e., (j) is cyclic), and can be parametrized in terms of 
the anisotropy parameters /3j obtained from expanding 
P^^ {oj, AE,9i) in Legendre polynomials P/(cos6'i), 

oo 

P^'iuj, AE, Oi) = P"'{uj, AE) Pj{oj, AE)Pj{cos 0i) . 

j=o 

(8) 

In the following, we will label the anisotropy parameters 
f3j, /3 = /32 and 7 = /34 (as e.g., in [18]). The integration 
over the angular part of one electron can be performed 
analytically. For two-photon double ionization from the 
ground state, the coefficients of Legendre polynomials 
with J > 4 vanish. In addition, odd multipoles /3j vanish 
because of parity conservation. Consequently, 

P^'{uj,AE,ei) = P^'{uj,AE)x 

X {l + /3(w,A£;)F2(cos6ii)+7(a;,A£;)P4(cos6li)} (9) 

where we have indicated the explicit dependence of the 
angular distribution and anisotropy parameters on the 
photon energy huj. 

For uncorrelated sequential emission with the ground 
ionic state as intermediate state, i.e., if each electron in- 
dependently absorbs one photon from a Is state, we ex- 
pect a dipole-like cos^(0) distribution for both electrons. 
This corresponds to /3 = 2 and 7 = 0. Deviations from 
this scenario, in particular correlated and temporarily 
confined joint emission manifests itself by deviations from 
these values. 

In one-photon DI of He, for example, the photon can 
only be absorbed by one electron, while the second elec- 
tron is mainly released due to shake-off. Experiments 
(cf. [46] and references therein) showed that indeed /3 is 
approximately zero for the slow electron, showing a clear 
sign of the isotropic shake-off process. Note that 7 is 
always zero for one-photon processes. For fku close to 
the one-photon DI threshold {huj > 79 eV), f3 is theo- 
retically predicted to be close to /3 « — 1 [47[ and also 
experimentally found to be negative (see e.g., [48[ and 
references therein), indicating a preference for emission 
perpendicular to the polarization axis. 

For two-photon DI a completely different scenario pre- 
vails. For photon energies well below the threshold for 
sequential ionization, the energy distribution is flat, sim- 
ilar to that of the one-photon process. However, the 
angular distribution is completely different. This is to 
be expected because of the different number of photons 
absorbed, and, thus the different amount of angular mo- 
mentum transferred. The one-electron angular distribu- 
tion at equal energy sharing, P^^ {AE = 0, 6*1), for differ- 
ent photon energies (Fig. 2b) closely resembles a Hertz- 
dipole distribution (~ cos^ 6, thin black line) suggesting 
that each electron absorbs one photon. Although the 
electrons have to strongly exchange energy to reach equal 
energy sharing, this interaction does not leave a visible 
trace in the one-electron angular distribution. The one- 
electron angular distribution does not provide detailed 



information on the correlated emission, unlike observ- 
ables related to the joint two-electron distribution, such 
as the triply differential cross section (TDCS) (see [21, 30[ 
and references therein), the forward-backward asymme- 
try between the electrons [22[, or the angle 6*12 between 
the electrons (as will be discussed below, see Fig. 6). 

The energy-differential one-electron angular distribu- 
tion characterized by /3(w, AE) and 7(0;, AE) for differ- 
ent photon energies huj between 42 eV and 80 eV (Fig. 2a) 
displays a remarkably similar dependence on AE but, 
unlike the energy distribution, also noticeable variations 
with huj. 

Deviations from the uncorrelated limit (/3 = 2, 7 = 0) 
are most pronounced for slow electrons {huj — 42 eV, with 
only 2.5 eV emission energy per electron above the double 
ionization threshold). Hence, emission perpendicular to 
the polarization axis can be observed with 7 = 0.6 and 
/3 = 1.6. As mentioned above, the anisotropy parameters 
are not symmetric relative to AE = 0. The reason is 
that we observe the angular distribution of the "first" 
electron, P°\Ei - £^2, 6*1). For AE < 0, we observe the 
slower electron, while for AE > the faster one. The 
electron repulsion force, i.e., the acceleration due to the 
interaction with the other electron, is the same for both 
electrons. However, the relative momentum change due 
to the interaction is much larger for the electron with the 
smaller momentum. 

For somewhat higher photon energies {hcu > 48 eV) 
the distribution mainly approaches the Hertz-dipole like 
cos^ 9 distribution with little difference between the 
NSDI and SDI regime highlighting the continuity across 
the SDI threshold also in the angular distribution. Fur- 
thermore, the distribution becomes approximately sym- 
metric {AE —AE) for larger huj. In the sequential 
regime, an asymmetry in AE remains visible but con- 
fined to large asymmetric energy sharings (|Ai?| — >■ Etot) 
where one electron is again slow. However, the ionization 
yields become very small for this case (see Fig. 1). 

We focus in the following on the anisotropy param- 
eters for mean photon energies above the DI threshold 
{huj > 79 eV) and consider both "long" pulses (T = 4.5 fs) 
and attosecond pulses (T < 500 as). For pulses of a 
few femtoseconds duration dynamically induced Fano 
resonances due to interferences between SDI via the 
ionic ground state and shake-up intermediate states are 
present [12, 22[. Note that in the limit of both ultrashort 
and long pulses these interferences disappear as one of 
the two pathways becomes dominant (the nonsequential 
background for T — > 0, the sequential shake-up peak for 
T^oo). 

The Fano resonance structure is also visible in the 
one-electron angular distribution. The two-dimensional 
one-electron energy- and angle-differential distribution 
P^^ {AE ,0i) (Fig. 3) reveals the resonant variation with 
AE at fixed 9i. At 9i — n /2, the distribution has a nar- 
row local maximum at energies where the sequential ion- 
ization process with shake-up in the intermediate state 
contributes. Outside the shake-up resonance position. 
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FIG. 3. Two-dimensional energy-angular differential distribu- 
tion P°'{AE, Oi) for an XUV pulse with T = 4.5 fs and hoj = 
80 eV as a function of the energy difference AE = Ei — E2 
between the two electrons, and emission angle 9\, relative to 
the polarization axis. The vertical dashed white lines show 
the expected positions of the peaks for the sequential process, 
see Fig. 1. The horizontal nodal line 9\ = 7r/2 is visible except 
near shake-up resonances. 



the persistence of the nodal Hnes of the Hertz dipole at 
9i = 7r/2 (see Fig. 2b) is clearly visible. 

Significant deviations can be only observed for energy 
sharings that correspond to resonant shake-up interme- 
diate states (n = 2, n = 3) for which higher angular mo- 
mentum (I > 0) states can be populated. Thus, strong 
angular momentum mixing and deviations from an s-like 
initial (or in this case, intermediate) state can occur. In 
turn, outside these shake-up resonances, the conserva- 
tion of the nodal line of the Hertz dipole at 61 = 90° is 
pronounced. This observation suggests that the /3 and 
7 parameters as a function of AE (or Ei, used in the 
following) are strongly interrelated. Enforcing vanishing 
emission probability at 6 = Tr/2 implies 



7'(i?i) = 



8 fl3{Ei 



3 V 2 



- 1 



(10) 



The symmetry-enforced value of 7'(i?i) using the cal- 
culated values of P{Ei) as input according to Eq. 10 
agrees with the directly calculated 7(i?i) remarkably 
well, in particular in between the n — 1 sequential main 
peaks (Fig. 4a). Exactly at the sequential peak we ob- 
tain j3 — 2 and 7 = as expected from the independent- 
particle process where the angular distribution is an un- 
correlated product of cos^ 9. Outside the main sequen- 
tial peaks with increasingly asymmetric energy sharing 
the angular distribution becomes elongated along the po- 
larization axis which is reflected in larger discrepancies 
between 7'(i?i) and ^{Ei). The elongated emission pat- 
tern can be qualitatively explained by the fact that highly 
asymmetric energy sharing can be reached through post- 
collision interaction (cf. [49, 50]), which is most effec- 
tive when both electrons are emitted coUinearly along 
the laser polarization axis. The high-energy electrons 



are thus expected to belong to the subset of electrons at 
small angles, explaining the elongated angular distribu- 
tion. 

For extremely unequal energy sharing, we find a strik- 
ing change of the anisotropy parameters at the positions 
where sequential ionization through an intermediate 
shake-up state is possible (Fig. 4e,f). These peaks and 
dips can be observed up to the n = 4 shake-up with the 
current pulse parameters. They are due to the fact that 
the intermediate ionic state for the sequential process 
can have higher angular momentum (Z 7^ 0). Hence for, 
e.g., shake-up to np states, the outgoing two-electron 
wave is not even approximately described by a (p,p) 
wave, but is dominated by the (s, d) contribution. 
The interference between different intermediate state 
channels leading to the same final state gives rise to 
the complex angular distribution (Fig. 4 (e) and (f) for 
the shake-up to intermediate states n = 2 and 11 = 3, 
respectively) . 

The present calculation of anisotropy parameters for 
TPDI can be compared with two previously published 
results for long pulses [19] and for attosecond pulses [18]. 
Comparison between the long-pulse limit of our present 
calculations (converged results are reached for T = 4.5 fs) 
and recent time-independent results by Ivanov et al. ]19], 
(Fig. 5a), show reasonably good agreement near equal 
energy sharing. However, because of the limitations and 
approximations in the convergent close-coupling (CCC) 
method ]19], such as the inequivalent treatment of the 
slow and fast electron, the intricate structures due to 
shake-up interferences appear to be missing. 

For ultrashort attosecond pulses dynamical resonance 
structures are washed out. At the same time, the strong 
temporal confinement of the two electron emission events 
induces strong electron correlations and renders the dis- 
tinction between sequential and nonsequential DI largely 
obsolete. The measure of the temporal confinement to 
within a few hundred attoseconds on the angular distri- 
bution is therefore of particular interest. We compare 
the one-electron angular anisotropy parameters for an 
XUV pulse with fiuj = 91.6 eV (previously investigated 
in ]18, 20, 23]) and a pulse duration of T = 150 as and 
T = 450 as (FWHM of the Gaussian pulse envelope) with 
previous results by Barna et al. [18]. At the sequential 
peak we find rather good agreement (Fig. 5b). However, 
for equal energy sharing the calculation by Barna et al. 
overestimates correlation effects. This is probably due to 
two reasons: (i) the angular distribution in ]18] was de- 
termined right at the end of the pulse when the electron 
momenta have not yet converged to their asymptotic val- 
ues, and (ii) only single electron angular momenta up to 
h,ma.x = ^2, max = 2 whcrc included, which is not suffi- 
cient for converged angular distributions (see [30] for a 
detailed study of the convergence of TDCS with angular 
momentum) . The most striking difference is the absence 
of the nodal line at 9 = Tr/2 in ]18] which we find to be 
well preserved for attosecond pulses. The preservation 
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FIG. 4. (a) Anisotropy parameters P{Ei) and 7(-Ei) as 
well as 7'(i?i) (dashed line), which according to Eq. 10 
enforces vanishing probability of electron ejection at 90° 
to the polarization axis, (b)-(f) show angle-resolved one- 
electron distributions at different electron energies, nor- 
malized to a value of 1 for 9 = 0°: (b) equal energy shar- 
ing, (c) asymmetric energy sharing, (d) main sequential 
peak (intermediate state: |He^ls)), (e) second sequential 
peak (intermediate state: |He^2Z)) (f) third sequential 
peak (intermediate state: |He"'"3Z)). Parameters: central 
frequency Hlu = 80 eV, duration T = 4.5 fs. The thin 
black line in (b)-(f) shows a cos^ distribution. 
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FIG. 5. Asymmetry parameters /3 (upper group of lines and squares) and 7 (lower group of lines and circles) for TPDI in the 
sequential regime with an XUV pulse with (a) huj = 90 eV in in the long pulse limit, T = 4.5fs (FWHM - sin^ pulse envelope) 
in comparison with recent time-independent results by Ivanov et al. [19] and (b) hui = 91.6 eV for ultrashort pulse durations, 
T = 150 as and T = 450 as (FWHM - Gaussian pulse envelope) in comparison with the calculations of Barna et al. [18]. 



of the nodal line suggests that for the limit of ultrashort gree of approximation, intact. 

pulses, TPDI proceeds by the absorption of one photon 

per electron as we have previously observed [21]. Strong 

interaction between two almost simultaneously outgoing 

electrons leads to strong energy redistribution but leaves 

the nodal line in the angular distribution, to a good de- 
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FIG. 6. Angular correlation characterized by the expectation value of the angle between the emitted electrons, {cos6i2){AE), 
for photon energies between 42 eV and 80 eV (the other field parameters are the same as in Fig. 1 and 2). 



V. ANGULAR CORRELATION 

Correlation effects in the spatial anisotropy of TPDI 
are expected to be more pronounced in two-electron 
observables derived from the joint probability density 
P^^ {Ei,rti, E2,i^2)- Following the observation that the 
energy spacing AE governs the electronic dynamics both 
below and above the SDI threshold (Eq. 5, Eq. 7) we first 
integrate over -Etotai = Ei + E2 to arrive at the reduced 
probability P^^ {AE, ili, f22) as a function of the energy 
sharing AE. As a measure of the angular correlation 
between the two outgoing electrons we choose to investi- 
gate the angle 6*12 between the electrons. The expecta- 
tion value (cos^?i2)(A£') provides a suitable quantity to 
measure the asymmetry in the joint angular distribution. 
A negative value of (cos 612) indicates back-to-back emis- 
sion, while for positive (cos0i2) the electrons are emitted 
into the same hemisphere. 

This measure displays, as a function of AE, very 
similar features over a wide range of photon energies, 
42 eV - 80 eV, for a pulse duration of T = 4.5 fs (Fig. 6) . 
Regardless of the photon energy, the electrons are 
strictly emitted back-to-back between the main peaks of 
sequential ionization with the ionic intermediate ground 
state He'^(ls) and in the same direction for energy 
sharings outside the main peaks at ±(Ji — 12). For 
photon energies below the SDI threshold (cos0i2)(Ai?) 
is strictly negative implying that electrons are pre- 
dominantly emitted back-to-back in the NSDI regime. 
Despite the apparent similarity of the AE dependence of 
the asymmetry (Fig. 6) and the anisotropy parameters 
(Fig. 2, Fig. 4) we emphasize one key difference: the 
behavior of (cos 0i2)(Ai?) is controlled by two-particle 
correlations while I3{AE) and 'y{AE) are one-electron 
variables. Unlike the anisotropy parameters, in the 
strictly sequential regime T — 00 when the emission 
is confined to the peaks of on-shell ionic intermediate 



states, (cos 6*12) (A£') would approach zero. The onset of 
this uncorrelated limit can be observed for moderately 
long pulses (Fig. 6 and Fig. 7): at the energy difference 
AE corresponding to the sequential process with the 
intermediate ionic ground state the asymmetry goes to 
zero, (cos 6*12) (A£;) = 0. 

We conclude our discussion of the various distributions 
as a function of AE by comparing two different pulse du- 
rations in the long-pulse limit (T — 4.5 fs and T = 10 fs) 
for a photon energy huj = 80 eV. Fig. 7 shows the energy 
distribution divided by the pulse duration P^' {AE)/T, 
the anisotropy parameters /3, 7, and {cos 812) {A E) . All 
three distributions show that convergence with pulse du- 
ration is already reached for T — 4.5 fs over a wide range 
of energy sharings. This underscores the linear scaling 
of the energy distribution for nonsequential ionization, 
i.e., for all intermediate states that are off-shell. Only 
close to the peaks of the sequential process do the two 
distributions deviate from each other. At the position of 
the resonances the yield scales quadratically with pulse 
duration, leading to a growth of the sequential contribu- 
tion with respect to the nonsequential background. An 
interesting consequence of this can be observed for the 
n = 3 shake-up interference in Fig. 7 (see inset) where 
the dip in the Fano resonance for T = 4.5 fs turns into a 
peak for T = 10 fs. 

Similarly, the (joint) angular distributions are con- 
verged in the purely nonsequential regions. In contrast, 
the resonant features at the sequential peaks get more 
pronounced for the longer pulse since the sequential con- 
tributions scale as whereas the nonsequential back- 
ground scales linearly with T. For example, for a pulse 
duration of T = 10 fs, (cos 6*12) {AE) approaches the long 
pulse limit (cos 6*12) = at the shake-up sequential peaks 
for n — 2 and n = 3, while it is still non-zero for the 
shorter pulse. 
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FIG. 7. Energy distribution divided by the pulse duration 
P^' {AE)/T, anisotropy parameters P, 7, and (cos0i2)(Ai5) 
as a function of AE for two different pulse durations: T = 
4.5 fs (blue dashed line) and T = lOfs (red solid line) and 
hio = 80 eV. 



VI. TOTAL CROSS SECTIONS 

We finally address the consequences of the continuity 
across the SDI threshold for angle-integrated cross sec- 
tions for pulses with femtoseconds duration. Starting 
point is the shape function (Eq. 4, Eq. 5). The conti- 
nuity across the SDI threshold is explicitly incorporated 
by the assumption of constant (or smoothly varying) ma- 
trix elements (Eq. 5) . Accordingly, at a fixed value of T 
the shape function is a smoothly varying function of the 
final-state energy i^totai as well as of the photon energy 
fujj. This is the reason underlying the notion of the "vir- 
tual sequential ionization" appearing in the vicinity but 
below the SDI threshold [15, 16]. The only remaining 
manifestation of the SDI threshold in Eq. 4, Eq. 5 is the 
vanishing argument AEni = 0, AEf„ = in the complex 
exponential functions when huj reaches the ionization po- 
tential of the ionic ground state, i.e., the intermediate 
ionic ground state becomes on-shell. As fiuj approaches 
this value from below, P^^{AE) smoothly increases as 
AE — >■ ±(/2 — /i) and the energy mismatches AEfn — >■ 
and AEni — >■ decrease. The only apparent disconti- 
nuity between the nonsequential and sequential regimes 
appears in the limit T — > 00 as the contribution of the 



shape function around the on-shell transition scales as 
~ T^, while all off-shell intermediate contributions for 
hu both below and above the SDI threshold scale lin- 
early with T as T — >■ 00. The width of the region around 
AEfn — AEni — exhibiting quadratic scaling is pro- 
portional to 1 /T. In the nonsequential regime only parts 
of the shape function with AEfn, AEni 7^ are accessi- 
ble. The total yield P^^ = / P°^{AE)dAE thus scales 
linearly with T for long enough pulses (cf. [22]). 

Conventionally, one defines a generalized two-photon 
double ionization cross section as 



0-(2) 



lim „^ 



-P 



DI 



(11) 



with Teff.2 = 35r/128 for sin^ pulses and a two-photon 
transition. Eq. 11 assumes that P^' is approximately 
constant over the spectral bandwidth of the pulse. At 
the SDI threshold the asymptotic scaling of P^^ switches 
from ~ T to ^T'^ and P^^ thus depends strongly on the 
(mean) photon energy. The pulse should then be long 
enough to possess a sufficiently small spectral bandwidth 
to uncover the correct asymptotic scaling. By the same 
token, however, the variation of P^' as a function of hu! 
close to the threshold becomes more pronounced as the 
pulse duration becomes longer and the approach to the 
underlying divergence of P^^ [AE) for AE — )■ ±(J2 — /i) 
is increasingly uncovered. Therefore, the hallmark of the 
SDI threshold is a pronounced rise of P^^ and, in turn, 

of (7(2). 

An alternative way of extracting cross sections from 
a time-dependent calculation has recently been used by 
Palacios et al. [12]. While it removes the uncertainty in 
the total energy of the pulse, the energy resolution in the 
intermediate state is still determined by the pulse dura- 
tion. Close to threshold long pulses are thus still needed 
for accurate extraction of the cross section. In addition, 
it requires very precise energy resolution in P^^{Ei, i?2). 
For simplicity, we use Eq. 11 in the following. 

The total cross section for nonsequential TPDI 
strongly increases just below the SDI threshold, as antic- 
ipated above and observed in a number of recent studies 
]12, 15, 16, 30, 32]. This increase is a direct manifesta- 
tion of the continuity across the threshold and can be 
qualitatively explained as follows: because of the time- 
energy uncertainty, the system can transiently "borrow" 
the energy that is missing for the sequential pathway to 
be accessible, if the intermediate state (i.e., the state af- 
ter one-photon absorption) is only transiently occupied. 
In other words, while the final DI state is on-shell, the 
intermediate ionic state is off-shell but only barely so. 
The occupation time of the intermediate state scales with 
the inverse of the borrowed energy. For photon energies 
close to the SDI threshold, the system has to borrow 
only very little energy. In turn, the possible occupation 
time of the intermediate state corresponding to the se- 
quential process increases leading to a large total yield. 
This process has previously been called "virtual sequen- 
tial" TPDI ]15, 16]. Eventually, when the SDI threshold 
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FIG. 8. Total TPDI cross section obtained using sin'^ pulses 
with different durations T. Spectral distributions of the pulse 
are displayed for one data point for each T, showing the area 
over which the cross section obtained is effectively integrated. 
The intensity is 10^^ W/cm^ for all pulses. Angular momenta 
up to I/max = 3 for the total angular momentum and Zi^max = 
'2,max ~ 7 for the single angular momenta are included. The 
radial boxes had extensions up to rmax = 1400 a.u. for the 
10 fs pulses. 



photon sequential process, since we only take partial 
waves with L = and L = 2 into account when calculat- 
ing cr(-2) . To verify that excitation of highly excited ionic 
states is not wrongly interpreted as double ionization, 
we also checked the population in these states (which are 
occupied through two-photon single ionization). For the 
longest pulse used here (10 fs) and at a central frequency 
of htj = 53.5 eV, ions in the n = 1 state after single ion- 
ization can be resonantly excited into n = 7 and n — 8 
by the second photon, whereas the population of higher 
n states is strongly suppressed. We have checked on the 
stability of both the ionic Rydberg population as well as 
the double continuum population as a function of prop- 
agation time after conclusion of the pulse and prior to 
projection for up to 3 fs. The temporal stability indi- 
cates that an unintended inclusion of high-lying Rydberg 
states in the calculated double ionization probability can 
be ruled out. 

We therefore conclude that the rise of the general- 
ized two-photon cross section as the SDI threshold is 
approached is a physical consequence of the continuity 
across the threshold and signature of the onset of the 
"virtual" SDI process. 



VII. SUMMARY 



is crossed, both transitions can be on-shell and the inter- 
mediate state for the sequential process can be occupied 
for an infinitely long time. 

In Fig. 8 we show the total nonsequential TPDI cross 
section (Eq. 11), as extracted from pulses with a sin^ en- 
velope and durations T from 0.5 fs (comparable to the 
frequently used 10 cycle pulses in literature) up to 10 fs 
(corresponding to a total pulse duration of 20 fs). At 
some of the data points, we plot the spectral distribution 
of the laser pulses that were used to extract the cross sec- 
tion values. Obviously, the frequency dependence of (7(2) 
near the SDI threshold can only be mapped out when the 
spectral width does not overlap with the above-threshold 
region. Even for low photon energies the cross section is 
"blurred" when too short pulses are used. This leads to 
a wrong slope of the cross section, as can be seen by 
comparing the results for T = 0.5 fs with longer pulses. 
Fig. 8 also demonstrates unambiguously that the rise of 
the cross section close to the sequential threshold does not 
result from an unintended inclusion of sequential contri- 
butions to the total double ionization yield. The spectral 
bandwidth of the 10 fs pulse is small enough (w 0.3 eV 
FWHM) that virtually no contributions from the on-shell 
sequential process are present even at a central photon 
energy of 54eV. The spectral contribution |J^(a;)|^ with 
hu! > I2 is only 0.026% of the entire spectral content 
of the 10 fs pulse at 54 eV (2.39% for the 5.5 fs pulse). 
From the value of P^^ just above I2 and the residual 
weight |J^(a;)|^, the contribution to cr(2) at 54 eV can be 
estimated to be less than 0.42%. 

We can also rule out any contribution from the three- 



We have discussed consequences of the continuity 
across the threshold from nonsequential double ionization 
(NSDI) to sequential double ionization (SDI) at hw = l2- 
We have shown that the energy difference AE between 
the two outgoing electrons provides a suitable variable 
to display the common features of the NSDI and SDI 
processes on an energy- and angle differential level. 

In particular, the singly differential energy distribution 
as a function of the energy difference between the elec- 
trons, P^^{AE), agrees excellently over a wide range of 
energy sharings, irrespective of the photon energy. The 
shape of the energy distribution can be understood by 
a simple model based on second-order time-dependent 
perturbation theory. 

We have also presented, for the first time, fully con- 
verged anisotropy parameters for TPDI over a wide range 
of photon energies (42 — 80 eV) and of pulse durations 
(150— 10000 as). An approximately conserved nodal line 
at 90° relative to the polarization axis, which remains 
undistorted by final-state electron-electron interaction, 
points to the fact that the ionization process is domi- 
nated by the absorption of one photon by each electron. 
Deviations from this emission pattern could only be ob- 
served in regions where one of the electrons is so slow that 
its direction can be easily changed, i.e., for very unequal 
energy sharing outside the main sequential peaks, or for 
low photon energies only slightly above the threshold for 
TPDI. A second notable exception is the region of dy- 
namical Fano resonances where additional intermediate 
pathways via excited ionic states open which can have 
non-zero angular momentum. In such cases, the angular 
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distribution does not even approximately resemble that 
of a Hertz-dipole shape. 

Additionally, we have shown that a further conse- 
quence of the continuity across the threshold is a sharp 
rise of the energy- and angle integrated "generalized" 
cross section just below the threshold which can be 
viewed as the onset of the "virtual" sequential ionization 
channel consistent with the smooth approach to the on- 
shell intermediate state. 
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